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ABSTRACT 

The solar acoustic p-mode line profiles are asymmetric. Velocity spectra have 
c§ ! more power on the low-frequency sides, whereas intensity profiles show the op- 

posite sense of asymmetry. Numerical simulations of the upper convection zone 
^ , have resonant p-modes with the same asymmetries and asymmetry reversal as 

^ ! the observed modes. The temperature and velocity power spectra at optical 

depth Tcont = 1 have the opposite asymmetry as is observed for the intensity and 
velocity spectra. At a fixed geometrical depth, corresponding to < Tcont >= 1, 
however, the temperature and velocity spectra have the same asymmetry. This 
indicates that the asymmetry reversal is produced by radiative transfer effects 
and not by correlated noise. 



Subject headings: convection — methods: numerical — radiative transfer — Sun: 
oscillations 
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1. Introduction 

For almost a decade, it has been known that the power spectra of solar acoustic modes 
are asymmetric, velocity has more power on the low frequency side and intensity has more 
power on the high frequency side of the power maxima (e.g. Duvall et al. 1993). The 
asymmetry reversal between velocity and intensity is thought to be due to the correlated 
background noise contribution to the intensity power spectra (Nigam et al. 1998). It is 
debated whether the asymmetry reversal occurs in velocity (Roxburgh & Vorontsov 1997) 
or intensity (Nigam et al. 1998) or both (Kumar & Basu 1999). Theoretical models predict 
asymmetries that depend on the source depth and type (Kumar & Basu 1999; Georgobiani 
et al. 2000b). Roxburgh & Vorontsov (1997) considered the superposition of dipole and 
quadrupole sources; Nigam et al. (1998) used a combination of monopole and dipole terms; 
Kumar & Basu (1999) show that the asymmetry reversal could be triggered even by dipole 
or quadrupole sources alone. 

Simulations of the shallow upper layer of the solar convective zone have resonant acous- 
tic modes like the Sun. The emergent intensity and the velocity in the photosphere have 
asymmetric spectra with the opposite asymmetry (Georgobiani et al. 2000a). In this letter we 
calculate the temperature and velocity power spectra at the continuum optical depth t — 1 
and at the geometrical depth corresponding to < r >= 1. At unit continuum optical depth 
the velocity and temperature have opposite asymmetry, with the velocity having more low- 
frequency power and the intensity more high frequency power. At fixed geometrical depth, 
however, the velocity and temperature have the same asymmetry, more low frequency power, 
for the fundamental mode. These results indicate that the asymmetry reversal is caused by 
radiative transfer effects, and not by correlated noise. 

2. Model of the Solar Convection 

We use a three - dimensional hydrodynamic code of Stein & Nordlund (1998, and refer- 
ences therein) to simulate the upper layers of the solar convection zone. The computational 
domain covers 6 Mm by 6 Mm horizontally and 3 Mm in vertical direction, from 0.5 Mm 
above the r = 1 surface to 2.5 Mm beneath it. The model includes non-gray, LTE radiative 
transfer. Horizontal boundaries are periodic, while vertical ones are transmitting. The spa- 
tial resolution is 100 km horizontally and ~ 50 km vertically, with a finer grid interpolated 
for solving the radiative transfer equation. The radiation field is calculated by solving the 
Feautrier equation along a vertical and 4 straight inclined rays, after averaging the Planck 
function into four bins by wavelength sorted according to opacity (cf Nordlund 1982; Nord- 
lund & Stein 1990, 1991). Snapshots are saved at 30 s intervals. We have simulated 43 hours 
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of solar time. 



3. Calculation of Power Spectra and Phase Relations 

We compare power spectra of intensity, temperature and velocity in order to shed light 
on a cause of the p-mode line profile asymmetry and its reversal between velocity and in- 
tensity. Prom our simulations, we have coordinate- and time-dependent quantities: vertical 
velocity V{x, y, z, t), temperature T{x, y, z, t), emergent intensity I{x, y, t), etc. To obtain a 

power spectrum at a particular geometrical depth zq, we use V and T at z = zq. To obtain 
V and T at a particular optical depth, say, r = 1, we calculate r(x, y, z, t) and interpolate V 
and T to the height at which t{x, y, z.t) = 1 for each position {x, y) and time t. Wc choose 
Zo to be the depth for which < r >= 1. Clearly, V{x, y) or T{x, y) at a geometrical depth zq 
will not be the same as V{x, y) or T(a;, y) at t = 1: we see deeper in the cooler intergranular 
lanes compared to the hotter granules. We investigate if this makes a difference in the power 
spectra and somehow affects line asymmetries. 

Wc separate the oscillation modes into radial modes, for which we average our data 
horizontally, and nonradial modes, for which we multiply the data by the corresponding 
spatial sines or cosines and then average horizontally (for more details, see Georgobiani et 
al. 2000a). We then Fourier - transform these time strings to get power spectra and phase 
relations. In our simulations, the first nonradial mode (horizontal wavelength \h — Q Mm) 
corresponds to a harmonic degree I = 740, because the simulation box represents a small 
fraction of the solar surface. We investigate the behavior of the emergent intensity, plasma 
temperature and plasma velocity for this first non-radial mode of the simulation (see also 
Georgobiani et al. 2000a). 

4. Results 

The emergent continuum intensity (Fig. 1) is a good measure of the plasma temperature 
at local instantaneous continuum optical depth unity, T{x,y,t) = 1, as expected from the 
Eddington-Barbier relations. The temperature at r = 1 has the same spectrum (Fig. 2) and 
phase (Fig. 3) as the emergent intensity. These phase differences are essentially zero for all 
frequencies. The spectra of the emergent intensity and plasma temperature at r = 1 have 
opposite asymmetry to the plasma velocity (Fig. 4) as is observed (e.g. Duvall et al. 1993) 
and as was discussed by Georgobiani et al. (2000a). The spectrum of the temperature at 
fixed geometrical depth zq, corresponding to the average continuum optical depth unity, < 
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T >x,y,t— 1, (Fig- 5) is, however, rather different from its spectrum at local instantaneous unit 
optical depth (Fig. 2), with different asymmetries especially noticeable for the fundamental 
mode. The temperature at < r >= 1 has the same asymmetry as the velocity. 

5. Discussion 

What changes the asymmetry of the temperature spectrum between measuring it at 
local T — 1 and average < r >= 1? We analyze the first non-radial fundamental mode 
that has the most prominent asymmetry and agrees closely with the corresponding solar 
£ = 740 mode. Figure 6 shows the velocity and temperature profile of this mode at average 
continuum optical depth one. It is clear that the velocity and temperature have the same 
profiles. Figure 7 shows the velocity and temperature profiles measured at local optical 
depth one, which is where one would see them. The velocity spectrum is hardly changed, 
but the amplitude of the temperature fluctuations is reduced by more than an order of 
magnitude. The high temperature sensitivity of the H~ opacity obscures high temperature 
gas and alters the height at which the gas is observed. This reduces the magnitude of the 
observed temperature fluctuations, but this reduction is not as great on the high frequency 
side as on the low frequency side. Hence, the mode asymmetry is changed. 

Why is the reduction of the temperature fluctuations different at high and low frequen- 
cies? At the fixed geometrical height < r >= 1, the temperature fluctuations were larger 
on the low frequency side of the mode. This produces a larger opacity variation on the 
low frequency side of the mode (Fig. 8), which in turn leads to a larger variation in the 
height where local T{x,y,t) = 1 (Fig. 9). The radiation temperature we see is equal to the 
gas temperature at optical depth unity, according to the Eddington-Barbier relations. The 
phases of temperature and height of unit optical depth are such that where the temper- 
ature at flxed geometrical depth is largest we observe the temperature at greatest height 
(smallest z) (Fig 10). Since the gas temperature is decreasing outward, the larger variation 
in the location where the radiation originates on the low frequency side of the mode leads 
to a smaller temperature variation there, while the smaller variation in the location where 
radiation originates on the high frequency side of the mode leads to a larger temperature 
variation there. We therefore observe a larger intensity variation on the high frequency side 
of the mode than on the low frequency side of the mode and thus the asymmetry is reversed 
from that observed in the velocity (and temperature when measured at a fixed geometrical 
depth). 
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6. Conclusion 

We have found that the emergent intensity and the temperature spectra at local in- 
stantaneous optical depth unity have the opposite asymmetry to the velocity as is observed, 
while the temperature at fixed geometrical depth corresponding to average optical depth 
unity has the same asymmetry as the velocity. This indicates that radiation transfer plays a 
crucial role in the asymmetry reversal observed between the intensity and Doppler velocity, 
and that this reversal is not due solely to effects of correlated noise. The mode asymmetry 
has been shown to depend on the separation of observation location and source location for 
the case of a simple localized (^-function source (Kumar & Basu 1999; Georgobiani et al. 
2000b). T(r = 1) and T(< r >= 1) are observed at different heights. Oscillation induced 
opacity changes vary the location of radiation emission (r = 1) in a way that reduces the 
magnitude of the temperature fiuctuations and reverses their asymmetry. 

This work was supported in part by NASA grants NAG 5 9563 and NAG 5 8053, NSF 
grant AST 9819799, and by a grant of computer time from the National Center for Super- 
computing Applications (which is supported by NSF). Their support is greatly appreciated 
and was crucial to the completion of this work. 
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Fig. 1. — Spectrum of intensity, I^. The solid curve is boxcar smoothed over Au — 65/uHz. 
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Fig. 2. — Spectrum of temperature, T^,^ at r(a;, y, t) — 1. Again, the solid curve is smoothed. 
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Fig. 4. — Spectrum of velocity, V^, at t{x, y,t) — 1. Its spectrum is almost the same at fixed 
geometrical depth corresponding to < r >x,y,t— 1- 
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Fig. 5. — Spectrum of temperature , Tj,, at < r >= 1. The temperature at < r >= 1 has 
the same asymmetry as the velocity. 
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Fig. 6. — The temperature (dotted) and velocity (solid) spectrum for the first non-radial 
fundamental mode measured at < r >= 1. Here and after, diamonds represent unsmoothed 
temperature profiles, whereas stars correspond to the parameter on the left y-axis. The 
temperature and velocity profiles look very similar to each other. Curves are smoothed over 
32 yuHz. 
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Fig. 7. — The temperature (dotted) and velocity (solid) spectrum for the first non-radial 
fundamental mode measured at r = 1. The amplitude of the temperature fiuctuations is 
nonuniformly reduced across the mode peak. 



-14- 




Fig. 8. — The temperature (dotted) and opacity (solid) spectrum for the first non-radial 
fundamental mode measured at < r >= 1. The larger temperature fiuctuations on the low 
frequency side of the mode produce larger opacity variations. 
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Fig. 9. — The spectrum of the height of local optical depth one, Zr(x,y,t)=i, (solid) and 
temperature (dotted) measured at local r = 1. The location of local optical depth unity 
varies more on the low frequency side of the mode, where the opacity variation is larger 
due to the larger temperature fluctuations. This reduces the temperature fluctuations at 
local optical depth unity more on the low frequency side of the mode compared to the high 
frequency side and, in accordance with the Eddington-Barbier relations, leads to smaller 
intensity fluctuations on the low frequency side of the mode and larger intensity fluctuations 
on the high frequency side of the mode. 
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Fig. 10. — The phase (sohd hne) between the height of r = 1 and the temperature at 
< T >= 1. Phase PS 180° means that the height of unit optical depth is greatest (smallest 
z and lowest temperature) when the temperature at < r >= 1 is largest. Also shown is 
temperature at r = 1 (dashed). 



